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The one-dimensional optical polaron is treated on the basis of the perturbation 

theory in the weak coupling limit. A special matrix diagrammatic technique is 

developed. It is shown how to evaluate all terms of the perturbation theory for 

the ground-state energy of a polaron to any order by means of this technique. 

^ \ The ground-state energy is calculated up to the eighth order of the perturbation 

\ theory. The effective mass of an electron is obtained up to the sixth order of the 
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perturbation theory. The radius of convergence of the obtained series is estimated. 
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(N ■ I. INTRODUCTION 

o , 

Nowadays there is a stable interest in low- dimensional struct urestU~H. The one-dimensional po- 
laron problem is relevant in semiconductor physics, where with state-of-the-art nano-lithography 

O . gated organic polymer conductors^. The treatment of the polaron problem in quantum dots can 



it has become possible to confine electrons in one direction (quantum wires)Q and in linear conju- 



be found in Ref.i'i. 

Now there are a lot of theoretical works, where the polaron problem is investigated by means 
of the perturbation theory (PT)H3. The series of the perturbation theory is useful for verifying 
the approximate nonperturbative methods in the weak-coupling limit§H§. The PT for the N - 
dimensional polaron has been developed in Ref.§, where the perspectives of 1/iV-expansions are 
discussed. The technique of 1/iV-expansions has been developed later for the optical polaron 
in RefH Up to now the first three terms of the weak-coupling expansion for the ground-state 
energy of the bulk polaron have been calculated^ (see alsoE-Hll), as well as two terms for the 
surface polaron energy! and three terms of the wire polaron energyiEl. The investigation of the 
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convergence of the PT series of the bulk polaron can be found in Ref.H. 

In this paper the matrix diagrammatic technique is developed for an optical large polaron. 
This technique permits to evaluate any term of the PT, in principle. The ground-state energy of 
the one-dimensional polaron is calculated up to the eighth order of the PT by using the developed 
technique. The radius of convergence of the PT series is estimated by means of Cauchy-Hadamar's 
criterion with respect to the calculated terms of the series. 

In Section II the matrix diagrammatic technique is developed. In Section III the obtained 
results for the ground-state energy and the effective mass of an electron are given. The radius of 
convergence of the PT series is also estimated. 

II. THE MATRIX DIAGRAMMATIC TECHNIQUE FOR THE POLARON PROBLEM 

IN THE WEAK COUPLING LIMIT 



The Hamiltonian of the one- dimensional optical large polaron is given b; 

H = T + £<W fl * + E (X* < e ' lkX + Vh a k e tkx ), (1) 
z k k 

where u>k is a frequency of a phonon with momentum k, for the optical polaron u>k = ou does not 
depend on k and V4 = 2 1//4 (a;/L) 1//2 ; p, x are the momentum and space operators of an electron; 
ajj~, afc are the creation and annihilation operators of a phonon with momentum k; L is a normalized 
length; a acts as a coupling constant of the electron - phonon interaction. Our units are such that 
h, uj and the electron mass are unity. Below we shall make usual simplifying assumption, such 
that the crystal lattice acts like a dielectric medium. It means that we can replace a sum J2k by 
an integral L J dk/2ir. 

Let us consider the weak coupling limit a C 1 for the polaron with the Hamiltonian ([!]). After 
doing Lee-Low-Pines transformation 

H' = U~ X HU, (2) 

m = u\&), (3) 

U = exp \i (P — ka^at) x], (4) 

k 

where P is a c-number representing the total system momentum, we obtain the Shrodinger equa- 
tion for Eq. (JTJ) in the next form 
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(H + Hi)\*') = E\*'), 

Ho = \{P -Y. ka k^kf + Y. a t a k, (5) 

/ k k 

Hi = Y. y k (4 + Ofc) • 
k 

Let us use the conventional perturbation theory^ for the SE Eq. (||) in the Fock basis, where H is 
an unperturbed Hamiltonian. Thus, if the zero-approximation of the vector state | ty') is a vacuum 
state of a phonon field |0) then the ground-state energy of H is given by Eq ^ = (0|H |0) = P 2 /2. 
It is easy to verify that all odd terms of the PT are equal to zero E^ = E^ = . . . = 0. The 
second order of the PT§ is 

Ef ] = {H x ±-H x ) 

2^ 2 a [™ difei 

(6) 



2tt J-ccE^ -i(P-^) 2 -l 
where ho = E^ — H and (0| . . . |0) — > (. . .), this term is defined by one connected diagram (see 
Eql.eps). The thick line corresponds to an electron propagation with momentum P — k\. The 
thin line corresponds to a propagation of a virtual phonon with momentum fci, a bold point on the 
thick electron line corresponds to a vertex, where one phonon either creates or annihilates. Below 
we shall give the Feynman rules for the connected diagrams represented in the matrix form. The 
quantity of the diagrams increases in the next orders of the PT: two connected diagrams in the 
forth order, ten connected diagrams in the sixth order, seventy-four connected diagrams in the 
eighth order and so on. There are also unconnected diagrams. These diagrams can be evaluated 
by differentiating the energy terms, which are the sums of the corresponding connected diagrams, 
with respect to E^ (see below). Note that all multidimensional integrals corresponding to the 
diagrams are evaluated by the residue theory. It can be seen from (|ToD . Thus, there is a technical 
problem to generate and evaluate all diagrams in the higher orders. 

Now we shall show how to build the matrix diagrammatic technique which permits us to 
generate all connected diagrams by means of any modern system of computer algebra (SCA). Any 
connected diagram can be represented in the n-th order of the PT by using n/2 by (n — 1) matrix 
||iV||, where n is an even number. For example, let us consider the energy term of the forth order 
(n = 4). It has the following forml! 

1 „\ 3 > 1 d r 1 - 2 



Ef = {H, ( i- H, ) ) + 



ho V ' 2 dE ^ 
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(fflT-ffl) 

ho 



(7) 



This term is defined by a sum of two connected and one unconnected graphical diagrams! (see 
Eq2.eps). These diagrams can be written in the following matrix form 



E, 



(4) 



1 1 




+ (1) 



0(1) 

oeP 



1 0, 

where z-th row of ||iV|| describes a history of propagation of i-th phonon with momentum h 
(i = 1,2. . .n/2). And j-th matrix column shows a distribution of phonons after passing j-th 
vertex, where one phonon either creates or annihilates (j = 1, 2 . . .n — 1). The value of = 1 
or corresponds to the existence or absence of i-th phonon between j-th and (j + l)-th vertex 
respectively. The generation of all connected diagrams for the n-th order is realized by selecting 
n/2 by n — 1 matrixes with respect to the next rules 

n/2 n _l 

A^ = 0orl, £i\^0, £%^0, 

i=l j=l 
n/2 n/2 



i=l i=l 
71-1 



(9) 



We have to keep only one arbitrary matrix among matrixes which are transformed into each other 
by permutating matrix rows. Thus, the whole set of connected diagrams can be got in the matrix 
form by means of any SCA. Using the graphical diagrammatic technique! it is easy to find the 
next rule for our matrix diagrammatic technique 



n-l 



(o) 



(2 l / 2 a\ n/2 r+™ r+°° 

— — / dh... dk n/2 
y 2n J J-oo J-oo 



E, 



1 



n/2 



n/2 



z i=l i=l 



(10) 



Any diagram represented in the matrix form corresponds to an analytical expression. So that in 
accordance with the rule (0) we have for (§) 



(1) 



2 l/2 



a 



1 1 1 

1 0, 

x 

1 1 

1 1 

x 



2tt 



(o) 



difei 
2a 2 



E^-UP-h 



-, -1 



dh 



E^-\{P-hf-l\ l 



E, 



-(P-h-k 2 ) 



(o) 



2a 2 
(2k) 



dh 



dh 



E 



E^-^P-h] 



(0) 



i -1 



E, 



-(P-h-k 2 ) 



-i 



E 



(0) 



- 2 (P - h) 



-1 



In order to summarize all unconnected diagrams in the n-th order we can use the general structure 
of the conventional perturbation theory series^. Note that the term of the PT 

n-l 

*«■ 

contains all connected diagrams. This term does not contain unconnected diagrams at all. The 



4 n)c = (Hx (^) ), (11) 



other terms of n-th order E^ n contain unconnected diagrams modifying the powers of the cor- 
responding electron propagators in the previous orders. If the dependence of Eq{s < n) on 
Eq ^ will be conserved explicitly then Eq can be represented as a function of Eq and their 
derivatives. For example, Eq 1 ^ is written for some particular follows 

E ( ^ = E 2 E>, (12) 
E { o )n = ~(3») a E'{ + E 2 {E' 2 f + (E 2 E 4 )', (13) 
4 8)n = E 2 {E' 2 f + 3^(E 2 ) 2 E' 2 E'i + ^{E 2 fE' 2 " 

+E A (E' 2 ) 2 + 2—E 2 E 2 E A + (E 2 E 6 )' + E A E' 4 

+2E 2 E' 2 E' 4 + ±(E 2 fE'l, (14) 

where E n = Eq 1 ^ and a prime denotes a derivative with respect to Eq°\ All integrals ( [TOP are 
evaluated analytically by means of the residue theoryil without expanding them in powers of P. 
Then the effective mass of an electron is defined by the next formula 

1 d 2 E 



(15) 

p=o 



2m* dP 

We should like to note that the suggested matrix diagrammatic technique is acceptable for any 
iV-dimensional optical large polaron, but the rule (|i~Cf ) has to be generalized with respect to the 
Feynman rules for N- dimensional polaron!. 



III. RESULTS 

Let us carry out an asymptotic expansion of the ground-state energy up to the eighth order 
of the PT. At first, it is necessary to generate and evaluate all connected diagrams for the corre- 
sponding orders with respect to the conditions and rule flTU|). At second, we have to summarize 
obtained diagrammatic terms and unconnected diagrammatic terms defined by Eq. (O)-(n). 
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Thus, the polaron ground-state energy up to the eighth order is E (P) = E^o^o""^), where 
the energy terms are defined by 



E<°\P) 

4 2 \p) 

E$\p) 



p2 

2 ' 



2 l/2 



(2 - P 2 ) 



— a = -a - — a + o(P ), 



(16) 
(17) 



P 2 (P 2 -4) + 6 



P 2 (P 2 -3) +4 



(2-P 2 )3/2( 4 _p2)l/2 (2-P 2 )' 

'3V2 



a 



1 a 2 + 



P 2 
32 



-5V2) a 2 + o(P 4 ), 



P< 6) (P) 



Sv^ 16 V 3 I 



P< 8) (P) 



1 = H \ 11472V2 a 3 + o(P 4 ), 

2^ 4 32v/2 96 V 6 / 

/ 442369 218861 151925 261335 \ , ,- s 
1 1 a 4 + o(P )• 



(19) 
(20) 



V 15456 7728^ 2208^ 2208\/6 / 
Since the terms Eq(P) and Eq 8 \p) are too bulky, we have only written out their expansion in 
powers of momentum P. Using Eq. (pf)-(l2"0D the ground-state energy of a slow-moving polaron 
is written as 

P 2 



Eo(P) 



2m* 



- a - 0.06066017 a' 



0.00844437a 3 - 0.00151488 a 4 + o(a 5 ). 



(21) 



The effective mass of an electron is defined by equation Eq. (O 



, a 5-V2 . 



33 183 



8 32^ 



+ 



1 98593 23<P\ 3 . 4 . 1 n r 
-J = a 3 + o(a 4 ) ~ 1 + 0.5 a 

2 V 13824 24^/ 



+0. 1919417a 2 + 0.0691096 a 3 + o(a 4 



Now let us compare the obtained asymptotic formula for the polaron ground-state energy with 
the energy obtained in the frame of the Feynman polaron theory0@ (see Table |). For a < 3.4 
the asymptotic energy Eq. (|21|) lies lower than the Feynman variational result Eq with maximum 
deviation about 4 percent. For a ^ 5 formula Eq. (^) is not correct because of the radius of 
convergence of the series is R ~ 5 (see below). Note that the first three terms of the energy Eq. 
(PTI) coincide with the same terms from Ref.!0. 



Let us estimate the radius of convergence of the PT series for E (0). The radius of convergence 
R can be evaluated by Cauchy-Hadamar's criterion^ 

R = lim R n = lim f lE^l/a^Y^ . (23) 

It is clear from Table [TT] that there is quite fast convergence of the sequence {R n } near the point 
a ~ 5. So if the nonevaluated higher order energy terms conserve the existent tendency to 
convergence of the sequence {R n } then the series Eq. has a finite radius of convergence 

R ~ 5. 

IV. CONCLUSION 

The main purpose of the paper is to develop the matrix diagrammatic technique for the optical 
large polaron problem in the weak coupling limit. The first four terms of the ground-state energy 
and the first three terms of the effective mass of the one- dimensional polaron are evaluated by 
means of this technique. The suggested technique is acceptable for any iV-dimensional optical 
large polaron. The obtained results are compared with the results from the Feynman polaron 
theory. The radius of convergence of the PT series for the one-dimensional polaron is estimated 
by Cauchy-Hadamar's criterion. 
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TABLE I. The ground-state energy Eq(P). 
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a 



-Eq(0) from (P) 



0.1 




0.100376 




0.100615 


0.5 




0.510063 




0.516315 


1.0 




1.044445 




1.070619 


1.5 




1.613146 




1.672654 


2.0 




2.236957 




2.334434 


2.5 




2.959682 




3.070245 


3.0 




3.828595 




3.896646 


3.3 




4.426768 




4.443709 


3.4 




4.639049 




4.635570 


3.5 




4.857770 




4.832468 


4.0 




6.047798 




5.898815 


4.5 




7.398112 




7.119062 


5.0 




8.908301 




8.518858 


TABLE II. First four terms of the sequence {R n } 


n 


2 


4 


6 8 


oo 


Rn 


1 


4.060207 


4.910708 5.068795 


R 
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